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The production of mesons in Au+Au collisions at RHIC and their propagation in a hot 
and dense nuclear medium is studied within the microscopic quark-gluon string model. The 
inverse slope parameter of the transverse mass distribution agrees well with that extracted from 
the STAR data, while the absolute yield of is underestimated by a factor 2. It appears that 
the fusion of strings alone cannot increase the yield either. Less than 30% of detectable 0's 
experience elastic scattering, this rate is insufficient for the full thermalization of 0. The directed 
flow of <j) at \y\ < 2 demonstrates strong antiflow behavior, whereas its elliptic flow rises up to 
about 3.5% in the same rapidity interval. As a function of transverse momentum it rises linearly 
with increasing pt, in agreement with the STAR data, and saturates at pt > 2 GeV/c. 

The investigation of nuclear matter under extreme temperatures and densities, and 
^ ■ the search for a predicted transition to a deconfined phase of quarks and gluons, the 
so-called Quark-Gluon Plasma (QGP), is one of the main goals of heavy ion experiments 
at ultrarelativistic energies. Both theorists and experimentalists are looking for genuine 
QGP fingerprints, that cannot be masked or washed out by processes on a hadronic level. 
The meson was proposed about two decades ago |I],|2| as one of the most promising 
QGP messengers because of the following reasons: Firstly, an enhancement of ss pairs 
in the QGP phase should lead to an enhancement of mesons, as well as other strange 
hadrons, while the production of in hh interactions is suppressed due to the OZI rule [Q. 
Secondly, since the interaction cross section with non-strange hadrons is small, will 
keep information about the early hot and dense phase. Thirdly, the meson spectrum 
is not distorted by feeddown from resonance decays. Finally, among other channels 
decays into a pair of kaons and, more rarely, into a lepton pair. Both channels have been 
explored for the reconstruction of the yield in Pb+Pb collisions at SPS energy. It was 
found that the multiplicity extracted from the di-kaon data || is significantly smaller 
than that obtained from the dimuons |4| (the so-called meson puzzle). A similar effect 
has been observed recently by the PHENIX Collaboration |5| in Au+Au collisions at 
RHIC. A possible solution of this puzzle is the scattering of at least one of the daughter 
kaons in the nuclear medium accompanied by in-medium modifications of kaon and 
masses. Note also that the width and mass of the meson are quite sensitive to the 
strangeness content of the surrounding medium. 
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In the present paper we study production, propagation, and freeze-out of mesons in 
Au+Au collisions at y/s = 130 AGeV within the microscopic quark-gluon string model 
(QGSM) [|K|. The transverse mass spectra of calculated for three different centrality 
bins are presented in Fig. [I] together with the experimental data reported by the STAR 

Collaboration pT| . Though the QGSM under- 
estimates the yield by a factor 2, it cor- 
rectly reproduces the slopes of all three dis- 
tributions. Results of the fit to an exponen- 
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Figure 1. m^-spectra of mesons 
in three centrality bins in Au+Au 
collisions at RHIC. 



Table 1 

Yields of 0, h~ , and ratio <f>/h~ at \y\ < 0.5 in 
Au+Au at RHIC calculated within the SFM. 



ted in Fig. |l| as well. Obtained values of the 
inverse slope parameter T are 378 (379+50) 
MeV, 349 (369+73) MeV, 370 (417+75) MeV 
for model (experimental) results, respectively. 
The absence of a centrality dependence in the 
m r spectra of may indicate the minor role of 
the final state interactions. The problem of 
a too low strangeness yield in the model can 
be solved by implementing the string fusion 
mechanism [O] . To check this, the yields of 0, 
h~, and their ratio 4>/h~ have been calculated 
in the string fusion model (SFM) Results 

are listed in Table [I]. Compared 
to the non-fusion case (NF), the 
string fusion (F) mode significantly 
decreases the and h~ multiplic- 
ities at midrapidity. An increase 
of the yield can be obtained by 
further rescattering (F+R) of sec- 
ondary hadrons via the channels 
KY -> 0iV, KE -> 0Y, KVt -> 0F, 
but it is still insufficient to match the 
data. The agreement between the calculations and the data can be improved by increasing 
the string tension or dropping the constituent quark masses H, by including the OZI- 
violation mechanism in the model (at low y/s this is the only channel for production), 
and by increasing the number of ss-pairs inside nucleons in multipomeron diagrams 0]. 

The mesons are produced either in primary nucleon-nucleon interactions via string 
fragmentation processes or in subsequent rescattering between produced strange hadrons. 
We found that in the first 4 fm/c the mesons are produced solely in string processes. The 
string channel dominates over the rescattering one up to 20 fm/c; in the later stages the 
rescattering channel prevails. This dynamical picture agrees well with that obtained for 
the production at SPS energy within the UrQMD model ||. The dN/dt distributions 
of 0, N, 7T, K, A, which are decoupled from the fireball after the last elastic or inelastic 
collision, are shown in Fig. ||. Compared to AGS and SPS energies [0, a substantial part 
of hadrons leave the system immediately after their production. Mesonic distributions 
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have sharp peaks at t = 8 — lOfm/c, while the distributions of baryons are wider due to 
the large number of rescatterings, which shift their dN/ dt maxima to later times. Figure 
D depicts the freeze-out distribution of the mesons over rapidity and emission time. 
We see that the bulk production of occurs in the central rapidity window \y\ < 1.5 
within the first 10 fm/c, while in \y\ > 2 areas the mesons are produced until the 
late stages of the system evolution: Because of the rapid expansion the central zone 
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Figure 2. dN/dt distribution of N,ir,K,A Figure 3. d 2 N/dydt distribution of the flnal- 
and over their last collision time t^ reeze . state mesons in the (y,t)-plane. 



becomes quickly dilute, its energy density drops, and early freeze-out of mesons, which 
have small interaction cross section with non-strange matter, takes place. According to 
QGSM, about 51% of detectable mesons decay after decoupling from the fireball. Less 
than 30% of mesons experience at least one elastic rescattering before the freeze-out. 
This relatively low rate of elastic collisions may not be sufficient for the thermalization of 
mesons, however, would it be enough for the development of noticeable anisotropic flow? 
- Other processes contributed to the flow evolution are production and absorption/decay 
of in dense medium. While the elastic collisions and resonance production of mesons 
increase the flow in the normal bounce-off direction, the absorption channel effectively 
reduces this normal component thus elongating the flow in the antiflow direction. 

Figure |] shows the rapidity dependence of the directed flow of 0, N, and K in minimum 
bias Au+Au collisions at RHIC. The slopes of the all distributions are negative at \y\ < 2, 
i.e. the antiflow component of the v\ dominates over its normal counterpart (see JH)| ' 



Similar antiflow slopes of vi(y) are developed by n and A [|T6| , p77[| ; its origin is traced to 



nuclear shadowing. At midrapidity \y\ < 0.5 the directed flow of is quite weak. 

The elliptic flow of mesons as a function of rj and p t is shown in Fig. [5] together with 



the model predictions for n and N , and recent experimental data |18| . The magnitude 
of the flow at midrapidity v\ (77) = 3.3 ± 1.1% is similar to that of pions and nucleons. 
At larger rapidity v% (y) is weaker than the flow of other hadrons. As a function of pt 
the elliptic flow of rises up to 2 GeV/c almost linearly with p t . The agreement with 
the experimental results is good, although large error bars do not allow more definite 
conclusions. The saturation and/or possible drop of the v^ipt) at higher p t is linked to 
the lack of rescattering in spatially anisotropic finite nuclear matter. 

In summary, the QGSM successfully describes m r spectra of in different centrality 
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Figure 4. Directed flow vi(y,all pt) of 
</>, iV, K in minimum bias Au+Au events at 
= 130 AGeV. 



Figure 5. The same as Fig. § but for the 
elliptic flow V2(y,all pt) and V2(pt,\v\ < 1) 
of 0, iV, if. Data are taken from MM . 



bins in Au+Au collisions at yfs = 130 AGeV, but underestimates the total <fi production. 
This problem cannot be cured alone by the fusion of strings. Mesons are frozen in the 
model earlier than baryons. The QGSM predicts unique antiflow slope at midrapidity for 
the directed flow of <fi, ir, N, K, A, and reproduces correctly the elliptic flow v^ipt)- 
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